Incoherent Soliton Excitations and Spin-Charge Separation in Blue Bronze 
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We present and analyze high resolution angle resolved photoemission (ARPES) data for K0.3M0O3 
(blue bronze). The observed Fermi surface, with two sheets that are one dimensional in nature, is 
consistent with a ladder material with only weak inter-ladder coupling. The marked broadening of 
the ARPES lineshape, a significant fraction of an eV, can be interpreted in terms of the intra-ladder 
interaction between the bonding and antibonding bands of the ladder that leads to spin-charge 
separation. The interaction also gives rise to asymmetry in the dispersions seen to the left and the 
right of a given Fermi point. With improved energy resolution, a high energy feature is revealed for 
the first time in the photoemission spectra near the Fermi momentum. We believe this new feature 
is indicative of a higher energy bound state of soliton excitations on a ladder. 

PACS numbers: 74.81. Fa, 74.90.+n 



Systems of two coupled chains can be viewed as a first 
step in crossing over from one dimension, with all of 
its exotic physics of Luttinger liquid induced spin-charge 
separation, to higher dimensions. How these properties 
are transformed in the course of this crossover is why 
such systems (also called ladders) have been intensely 
studied by theorists. While theoretical progress on this 
problem has been considerable, with work in predicting 
rich excitation spectra, dynamical generation of spectral 
gaps, existence of preformed pairs (lTfl5l]. the progress on 
the experimental side has been slower both due to lim- 
ited numbers of materials with a suitable structure and 
limited experimental accuracy. 

K0.3M0O3 (blue bronze) is a quasi-one-dimensional 
(ID) material consisting of an array of weakly coupled 
pairs of conducting ladders [itjj . Band structure calcu- 
lations predict that the chemical potential is crossed by 
two slig htly warped bands corresponding to 3/ 4-filling 
[T7I Il8l |. These bands, arising from the uncoupled lad- 
ders, are composed of bonding (B) and anti-bonding 
(AB) bands with Fermi wave vectors Kfb and Kfab re- 
spectively. The Fermi surface nesting leads to charge den- 
sity wave (CDW) formation with a wave vector along the 
chains equal to Kfb + Kfab- A three-dimensional (3D) 
phase transition into a charge-ordered insulati ng p hase 
takes place at a temperature Tcdw = 180K (l9ll20l| even 
though the magnetic susceptibility starts to decrease well 
above the transition, experiencing a one-third drop in the 
interval between 700K and 180K [l9j. The correspond- 
ing activation gap is estimated to be 30-40 meV plj . 
Although the CDW related properties are well studied in 
blue bronze [22[, its characteristics as a ladder material 
are much less well understood. 

In this paper we report the results of high resolution 
angle-resolved photoemission (ARPES) measurements of 
the ladder material K0.3M0O3. With much improved in- 
strumental resolution, we resolve for the first time two 
features in the measured photoemission spectra. These 
two features are broad peaks in regions of the Brillouin 



zone near the Fermi vector, dispersing through energies 
of a significant fraction of an eV. We interpret them as 
ladder excitations originating from the same interaction 
that underpins the three dimensional charge order exist- 
ing in this material at low temperatures. The breadth 
of these peaks arises from the presence of gaplcss charge 
excitations (holons) and gapful spin solitons (spinons) 
with markedly different velocities. Assuming that the 
electron-lattice interaction plays a secondary role and the 
physics is electron-driven, we obtain a good description 
of the higher energy portion of the spectral function with 
two SU(2) Thirring models. 

The K0.3M0O3 single crystals were grown by an elec- 
trolytic reduction method [23| . High resolution angle- 
resolved photoemission measurements were carried out 
with a Scienta R4000 electron energy analyzer 24|. For 
the band structure measurements and Fermi surface 
mapping (Fig. 1), we used a helium discharge lamp as the 
light source with a photon energy of hv =21.218 eV. The 
overall energy resolution was set at 10 meV and the an- 
gular resolution was 0.3°, corresponding to a momentum 
resolution of 0.0091 A -1 at the photon energy of 21.218 
cV. For the high-precision ARPES measurements (Fig. 
2 and Fig. 3), a vacuum ultra-violet (VUV) laser with a 
photon energy /w=6.994 eV was used as a light source. 
The energy resolution used in this case was 1 meV. The 
angular resolution is 0.3°, corresponding to a momentum 
resolution of 0.004 A -1 at the photon energy of 6.994 
cV. The Fermi level is referenced from measuring a clean 
polycrystallinc gold that is electrically connected to the 
sample. The crystal was cleaved in situ and measured in 
vacuum with a base pressure better than 5xl0 -11 Torr. 

Fig. la shows the bare Fermi surface (FS). The data 
for this figure are collected at 80K. The FS is open and 
consists of weakly warped sheets, indicating a weak inter- 
ladder coupling. Based on the the warping of the ex- 
tracted Fermi momentum (blue dots in Fig. la), one can 
estimate the coupling to be t± sa 100 meV . Fig lb shows 
the band structure measured along the T-Y cut at 80 K. 
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FIG. 1: a) The bare Fermi surface of K0.3M0O3 extracted 
from the measurements. Data are symmetrized with respect 
to the k||=0 line. Blue and red dots stand for the extracted 
Fermi momentum along the Fermi surface sheets. The CDW 
vectors are marked by the dashed arrows, b) The measured 
band structure along T-Y (kd = 0) showing two bands, the 
bonding band (B) and the anti-bonding band (AB), crossing 
the Fermi level. The solid lines mark bands predicted from 
band structure calculations fl8l |. 



At energies less than leV we can clearly distinguish two 
bands (B and AB). The non-interacting bands cross the 
chemical potential at 0.97 7r/b, 0.6 7r/b, -0.55 7r/b and - 
0.96 7r/b, consistent with the band structure calculations 
111 Ref. hj and previous ARPES measurements [25J, [26j . 
The Fermi surface nesting vector is thus estimated to be 
w 1.52 7r/b, also in good agreement with the measured 
CDW vector [13, S^. 

In Fig. 2 and Fig. 3, we present high resolution laser 
ARPES data for the band structure along the ladder at 
k c [ = 0. Due to the effects of matrix elements, the B 
band is more clearly visible than the AB band (Fig. 2a). 
Its electronic structure exhibits two features clearly dis- 
tinguished in the second derivative image (Fig. 2b), seen 
as excitation branches near Kfb- The lower binding- 
energy branch was observed in earlier ARPES measure- 
ments [ID, , while the second higher energy branch has 
not hitherto been seen. The lower branch, starting at 
ss 60 meV below Ep, reaches in energy as low as -0.8 eV 
as fen varies from w/b to 0.4 ir/b (red circles in Fig. 2b). 
The second branch with a lower spectral weight starts 
at « 150 meV, remains distinct until roughly w —200 
meV, and then merges with the first branch (blue circles 
in Fig. 2b). The two branches gradually become indis- 
tinguishable when the temperature increases (Fig. 3(a - 
e)). At high temperature (240K), only one broad hump is 
observed (Fig. 3e). A further feature not previously seen 
but that we observe is the left-right asymmetry of the low 
energy branch: the slope of the band on the right side 
of the Fermi momentum (fen > Kfb) vr ~ 0.96 eVA is 
smaller than that on the left side (k\\ < Kfb)'- ~ 4.42 
eVA. These roughly correspond to the expected veloci- 
ties of the B and AB bands [18| ■ We will argue that both 
the presence of two excitation branches and the asym- 



metry of their dispersion can be explained in terms of 
one-dimensional intra-laddcr interactions. 

There are two standard approaches to investigate 
quasi-lD systems. One can take the transverse tunneling 
into account on the first step and treat the interactions 
via a perturbativc rcnormalization group and mean field. 
This approach predicts all sorts of Fermi surface insta- 
bilities and related 3D phase transitions (CDW, super- 
conductivity, etc.), but has difficulty in describing strong 
fluctuations and all related effects such as incoherence 
of the excitations [29j. An alternative approach is to 
consider the ID subsystem (in this ladder) with 

all of its interactions as a basis for calculations, use the 
powerful non-pcrturbativc techniques available in one di- 
mension, and then treat the intcrladder interactions as 
a perturbation [3014321 ] . The strong precursor effects ob- 
served such as the drop in the magnetic susceptibility [l9| 
and the absence of Fermi crossing in a wide temperature 
range above the CDW transition (Fig. 3) suggest that 
the second approach better captures the physics of the 
material in question. This point of view is further sup- 
ported by our ARPES measurements showing a strongly 
incoherent single particle spectral function (Fig. 2c). 

The previous attempts to explain the broad ARPES 
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FIG. 2: a) High resolution measurement of the band struc- 
ture along r~Y momentum cut with 1 meV energy resolution. 
Due to matrix element effects, only the B band is clearly ob- 
served, b) Second derivative with respect to energy of the 
data in (a). Two structures around ki?B, corresponding to in- 
coherent soliton and bound state (see text), are clearly seen, 
c) The corresponding photoemission spectra (EDCs). Two 
features are clearly resolved for the spectra near the Fermi 
momentum k.FB- Red (blue) circles mark in b) and c) mark 
the dispersion of the soliton and the bound state. 
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FIG. 3: a) - d) EDC second derivative images of the B band 
along r~Y momentum cut at different temperatures. The data 
are divided by the Fermi-Dirac function, e) Extracted Kfb 
EDCs at different temperatures. Two energy scales of inco- 
herent soliton and bound state are marked by blue and red 
dashed lines, respectively. 



spectrum have been based on small polaron theory 
[2ol [33|]. However, the shear size of the gaps and the 
wide energy range of the incoherence tails, spreading into 
a significant fraction of an eV, suggests that these fea- 
tures have an electronic origin. This incoherence cannot 
be explained in terms of the standard band structure 
theory even modified with 3D Fermi liquid corrections. 
However, there are non-perturbative calculations for the 
single electron spectral functions which we will argue pro- 
vide a qualitative fit to the data. 

In all ID theories the incoherence originates from quan- 
tum fluctuations of the gapless holon mode (in the case 
at hand it is the mode that describes the dynamics of the 
charge density wave). The spectrum of this mode at low 
energies has a linear dispersion, i.e., u = v c \k\. Since the 
CDW instability occurs at the wave vector kpB + kpAB, 
both bands are involved in creation of the order param- 
eter. However, as is clearly seen on Fig. 2b, these bands 
have very different Fermi velocities vfb and vfab - a 
feature that is not present in standard theories of 2-leg 
ladders where this difference is neglected. This velocity 
difference precludes the emergence of an SO(6) symme- 
try, the latter effect being one of the most interesting 
predictions of the low energy theory of ladders 0, Q- 
However we will argue that a remnant of this symmetry 
remains in the form of an electron- hole bound state. 

Below we describe a low energy theory which is chosen 
in such a way that (i) it maximizes the fluctuations lead- 
ing to the observed incoherency of the ARPES lineshapes 
and (ii) it gives rise to a zero temperature CDW instabil- 
ity at hps + kpAB which, in the presence of 3D coupling, 
is shifted to finite temperatures. Since the spectral gaps 
are much smaller than the bandwidth, it is reasonable to 
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FIG. 4: Schematic of interaction between pairs of Fermi 
points, (Kfb,—Kfab) and (Kfab,— Kfb)- For the first 
pair we show the asymmetry in the dispersion on either side 
of these two Fermi points. The width of the red stripe indi- 
cates the expected intensity of the dispersing incoherent soli- 
ton (gapped spinon and gapless holon) observed by ARPES. 
The blue stripe marks the corresponding intensity of a two 
particle continua composed in part by a dispersing bound 
state. The intensities at —Kfab are shown as lower because 
of matrix element effects. 



linearize the dispersion close to the Fermi points by intro- 
ducing slowly varying right- and left-moving components 
of the fcrmionic fields 
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°1> Lia) J=B,AB. (1) 



Because of the different Fermi velocities in the bonding 
and anti-bonding bands, we expect the theory will not 
have a symmetry higher then U(l)xSU(2). Moreover, 
we expect the strongest interaction to occur between the 
bonding and anti-bonding bands, foreshadowing the rise 
of the observed 3D CDW order. To this end we consider 
a model where the fermions at Kfb interact with the 
fermions at —Kfab and similarly those at —Kfb inter- 
act with fermions at Kfab- This leads us to describe 
the system as two decoupled ID theories, H = H\ + B-j- 
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with H2 = H\{AB -O- B) and where r a are Pauli matri- 
ces. The coupling constant g is equal to the Kfb + Kfab 
Fourier component of the rcnormalizcd Coulomb inter- 
action. We take it to be positive, corresponding to an 
attractive interaction that is necessary for CDW for- 
mation. The spectrum of both Hi and H2 consists of 
a gapless, spinless holon mode and a chargeless spinon 
mode with gap A ~ W(g/v) 1 ^ 2 cxp(— nxv/g") where v = 
{vfb + Vfab)/%- The other coupling go comes from for- 
ward scattering, the k = component of the Coulomb 
interaction. Both g and go are associated with a renor- 
malization of the spinon and holon velocities. 

The spectral function corresponding to this theory can 
be obtained from the one for the case of equal Fermi 
velocities v, studied in [30I ] by a simple Galilean transfor- 
mation: t' = t, x' = tAv + x, Av = (vfb — vfab)/^- 
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FIG. 5: EDCs for the B band for k- vectors ranging from k = 
0.47rb _1 to k = 1.37T& : a) before averaging for differential 
surface doping; b) after averaging with a Gaussian of width 
0.077Tb- 1 . 



We so obtain: 



(TfoMO, 0)^(0,0)) 



(e ikF ' x G Rj (r, -x) + e ik * 'C,., (r, x)) , 

j=B,AB 

G LA b{t,x) 



sj (v s+ t + ix)(v c+ r + ixj 



x exp 



— \/ (v c+ t + ix)(v c -T - ix) 



(3) 



where Z ~ 1, w c ,s.± = v c , s ± Aw, and Glb{t,x) = 
Glab{t,x,vab vb)- The corresponding retarded 
Green function in Fourier space is 



Glab (w, *) = F(w + *At>, fc) 

1 + a 

a) — 7J r fc 



F( W> ft) = 



A 2 + u, 2 £; 2 - a) 2 



A + v/A 2 + V 2 fc 2 -w 2 ) 2 



I- a 
1 + a 



(Q - v c ky 



-1/2 



(4) 



where a = v s /v c and Q = ui + id. 

The most notable prediction of this theory is that the 
spectral features at any Fermi point (±Kfb, ±Kfab) 
will appear to disperse with different velocities to its right 
and left. This is illustrated schematically in Fig. 4 for the 
pair of Fermi points Kfb and —Kfab- We also see the 
observed asymmetry on the spectrum in Fig. 5 where 
we plot the expected line shapes for a set of EDCs for 
different k's in the vicinity of Kfb'- for k < Kfb the 
spectral weight disperses with velocity vfb while for k > 
Kfb the spectral weight, while much reduced, disperses 
with velocity vfab- 

In Fig. 5a we see how spin-charge separation leads to 
broadened features in the ARPES lineshape giving rise 
to two singularities in the spectral weight at an upper 
and a lower energy that disperse with different velocities 
(the spin and charge velocities) as k moves away from a 



Fermi point. However, different spin and charge veloci- 
ties alone cannot explain the broad linewidths observed. 
We believe a part of the broadening of the lineshapes 
arises from K0.3M0O3 being a poor conductor. Poor con- 
ductivity will lead to ladders in different regions on the 
surface of the sample to be at different dopings. This 
in turn will lead to a spatial dependence of Kfb and 
Kfab- In fact, such dependence was observed in a scan- 
ning tunneling microscopy study of a sister blue bronze 
Rbo.3Mo03 [3J|, where Kp was found to vary across the 
surface from about 0.656* to about 0.806* [b* = 2n/b). 
This could also explain the significant O.lA^ 1 breadth 
of the momentum distribution curve at Ep in K0.3M0O3 
reported in Fig. 5(c) of [2^|. To account for this effect we 
convolve the lineshapes with a Gaussian weight function 
of width = 0.077T& -1 . We do this on the r.h.s. of Fig. 5 
where we observe a considerable increase in broadening. 
The ARPES measurements on the ID Mott-Hubbard in- 
sulator, SrCu02, considered the best evidence for spin- 
charge separation, also yield broad peaks much broader 
than the corresponding theoretical prediction (see Fig. 
3 in [1^). We thus believe that such broadening is an 
inherent feature of ARPES measurements on insulating 
systems originating, again, from some inhomogeneity or 
disorder present in the samples. 

We have so far ignored interactions between fermions 
living at Kfb,— Kfab and —Kfb,Kfab- However, 
such interactions are present even though they are weaker 
than those responsible for the 3D CDW order. In a 
ladder system with equal velocities in the bonding and 
anti-bonding bands, the theory governing the low energy 
response has a dynamically generated SO (6) symmetry, 
and so predicts that there will be a bound state of the 
spin degrees of freedom of the electrons 0] . This bound 
state contributes to a response above the gap, A, of the 
single particle spectral function. In a purely SO (6) theory 
this response leads to a feature in the spectral function at 
an energy u = (1 + \/2)A for k = Kfb, Kfab- While in 
a theory with broken SO (6) symmetry (i.e., vb ^ vab) 
we cannot predict where exactly this threshold will oc- 
cur, we can say that this bound state will survive even 
in the presence of symmetry breaking, a consequence of 
the bound state having significant binding energy mak- 
ing it stable to perturbations. This bound state offers an 
explanation for the second peak observed in Fig. 2b at 
an energy roughly twice that of the gap. 

In conclusion, we have presented high resolution 
ARPES data for the CDW material K .3MoO 3 . We have 
argued that the features in the ARPES lineshapes, in- 
cluding incoherent broadening and asymmetric disper- 
sions from the same Fermi point, can be understood pri- 
marily as arising from strong correlations in one dimen- 
sion. 
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